Neutrino Masses, Mixing and New Physics Effects 
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We introduce a parametrization of the effects of radiative corrections from new physics on the 
charged lepton and neutrino mass matrices, studying how several relevant quantities describing the 
pattern of neutrino masses and mixing are affected by these corrections. We find that the ratio 
uj = sin 0/ tan 9 at m is remarkably stable, even when relatively large corrections are added to the 
original mass matrices. It is also found that if the lightest neutrino has a mass around 0.3 eV, the 
pattern of masses and mixings is considerably more stable under perturbations than for a lighter or 
heavier spectrum. We explore the consequences of perturbations on some flavor relations given in 
the literature. In addition, for a quasi-degenerate neutrino spectrum it is shown that: (i) starting 
from a bi-maximal mixing scenario, the corrections to the mass matrices keep tan^atm very close 
to unity while they can lower tan#Q to its measured value; (ii) beginning from a scenario with a 
vanishing Dirac phase, corrections can induce a Dirac phase large enough to yield CP violation 
observable in neutrino oscillations. 

PACS numbers: 14.60.Pq,12.60.-i,11.30.Hv,11.30.Er 



I. INTRODUCTION 

Our present knowledge of neutrino masses and mixing 
is mainly provided by various neutrino oscillation experi- 
ments, which give us information on the two independent 
mass squared differences, as well as on the three angles 
characterizing the leptonic mixing matrix. In the future, 
the study of CP violation in neutrino oscillations may 
allow to determine the Dirac-type phase entering in the 
leptonic mixing matrix, and neutrinoless double p decay 
experiments may provide the value of the effective Majo- 
rana mass. Despite the great achievements of oscillation 
experiments, there is still much to be learned about neu- 
trinos. One of the most problematic issues in neutrino 
physics is the lack of information on the mass spectrum, 
since the mass squared differences do not fix the abso- 
lute scale of neutrino masses. Indeed, the spectrum can 
exhibit a strong hierarchy, as in the case of quarks and 
charged leptons, or on the contrary, be quasi-degenerate. 

At present, there are different extensions of the stan- 
dard model (SM) which propose mechanisms for the gen- 
eration of neutrino masses through the enlargement of 
the SM particle content. The addition of heavy right- 
handed neutrino singlets constitutes a simple and eco- 
nomical way to give left-handed neutrinos small masses 
through the seesaw mechanism 

Bll SB Another sim- 
ple possibility relies on the extension of the scalar sector 
with a heavy SU(2)r_ triplet, with or without supersym- 
metry (SUSY) Within supersymmetric mod- 

els, neutrino masses can also arise from i?-parity violating 
interactions [jj , where the atmospheric and solar neutrino 
mass scales are generated at the tree level and radiatively, 
respectively. Recently, a new supersymmetric source of 
neutrino masses and mixings has been found considering 
non-renormalizable lepton number violating interactions 
in the Kahler potential [lfj rather than in the effective 
superpotential. All the above scenarios predict the ex- 
istence of light Major ana massive neutrinos. Suppressed 



Dirac neutrino masses are unnatural in conventional the- 
ories since they usually require extremely small Yukawa 
couplings. However, this problem can be surmounted by 
endowing the particle content of the theory with extra 
fermion or scalar fields and/or introducing new symme- 
tries DH IH El El El Naturally small Dirac neutrino 
masses also arise in extra-dimensional theories as a con- 
sequence of the small overlap between the wavefunctions 
of the usual left-handed neutrinos in the brane and the 
sterile right-handed ones in the bulk (or in other branes) 

El El El 13 MM 

In addition to the mechanism for the generation of neu- 
trino masses, one of the most intriguing aspects of lep- 
tonic physics is the experimental evidence that two of the 
lepton mixing angles are large, in contrast to the small 
mixing observed in the quark sector. The deep under- 
standing of the neutrino mass suppression mechanism 
and the bi-large leptonic mixing constitutes one of the 
most challenging questions in particle physics. A theory 
of leptonic flavor should provide a plausible explanation 
for the bi-large mixing, as well as for the neutrino mass 
spectrum. At the same time, it should predict relations 
among these quantities. There have been in the literature 
a large number of suggestions in this direction, consisting 
in the introduction of flavor symmetries or the assump- 
tion of specific textures for the leptonic mass matrices 

Ei. 



The tree-level predictions of any theory of leptonic fla- 
vor are subject to higher order corrections, which are 
computable if the particle content and parameters of the 
new physics model are specified. Still, in the absence 
of a "standard theory" of lepton flavor, it is pertinent 
to ask ourselves about the possible consequences of the 
(unknown) radiative corrections on the tree-level predic- 
tions of the model. The aim of this paper is to investi- 
gate these effects, focusing on: (i) the modification of the 
tree-level values for the ratio of mass squared differences, 
mixing angles and CP violation parameters; (ii) the ef- 
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feet on various flavor relations. This task will be carried 
out in two steps. We will first propose a parametrization 
for the unknown radiative corrections to the lepton mass 
matrices, based on rather general arguments on weak- 
basis independence. In this parametrization, the formal 
structure of the corrections is fixed, up to undetermined 
complex coefficients. Then, the possible effects of the ra- 
diative corrections is investigated taking these coefficients 
as random, and performing a statistical analysis of the 
average behavior of the parameters and flavor relations 
under consideration. 

It is important to note here the difference between our 
framework and other studies in the literature regarding 
random neutrino matrices p3. I2M l2rl [27j ■ In these anal- 
yses, the neutrino mass matrices are taken as random at 
the tree level, with a subsequent discussion of the predic- 
tions for the neutrino mass spectrum and mixing angles. 
In contrast with this approach, in this paper we take the 
tree-level charged lepton and neutrino mass matrices as 
fixed by some theory of leptonic flavor, and consider ran- 
dom perturbations (from radiative corrections) to them 
[28l| . These corrections are not completely arbitrary: the 
different terms of the perturbations must have a deter- 
mined formal structure (which is dictated by weak basis 
independence) but have random coefficients. 

This paper is organized as follows. In Section [n] we 
briefly review the current status of neutrino oscillation 
data. The parametrization of the unknown corrections 
to the lepton mass matrices is derived in Section IlIII In 
Section HVI we present our results, discussing (i) the effect 
of perturbations on masses and mixing parameters; (ii) 
the influence on some flavor relations obtained for some 
specific patterns in the literature; (iii) the consequences 
of the corrections in some special limits of interest. Our 
conclusions are summarized in Section Ivl 



II. NEUTRINO OSCILLATIONS: PRESENT 
STATUS 

All presently available neutrino oscillation data can be 
accommodated within the framework of three mixed mas- 
sive neutrinos 29]. 1 The first KamLAND results |2jj se- 
lect the large mixing angle MSW solution as the only 
surviving explanation for the solar neutrino problem. In 
addition, dominant solar neutrino conversion based on 
non-oscillation solutions are now excluded [3ll l32l | . Re- 
cently, the Sudbury Neutrino Observatory (SNO) has re- 
leased the improved measurements of the salt enhanced 
phase [3i| which, together with all the solar and Kam- 
LAND neutrino data, allow for a better determination of 
the oscillation parameters. In particular, the high-Am 2 
region (Atoq > I0 -4 eV 2 ) is now only accepted at the 3a 



1 We will not consider here the results from the liquid scintillator 
neutrino detector (LSND). 



level and maximal solar mixing is ruled out by more than 
5<7, rejecting in this way t he p ossibili ty o f bi-maximal 
leptonic mixing IS El HI El E3 Concerning 
the atmospheric neutrino sector, the water Cherenkov 
Super-Kamiokande (SK) |4lj and long-baseline KEK- 
to-Kamioka (K2K) (43 experiments indicate that neu- 
trino flavor conversion due to neutrino oscillations in the 

— > v T channel provide by far the most acceptable and 
natural explanation for the observed disappearance. 

Regarding the absolute values of neutrino masses, the 
situation is not so satisfactory. At present, the most 
stringent direct bound on the neutrino mass is provided 
by the Mainz and Troitsk ^4| experiments, which 
have set a maximum value for m Ve of 2.2 eV. Neutrinoless 
double beta-decay measurements may also be valuable to 
disentangle the pattern of neutrino masses, although in 
this case certain subtleties have to be considered f4a| . 

Throughout this paper we adopt the neutrino mass 
ordering mi < < 7773 in such a way that, for the 
hierarchical (HI), inverted- hierarchical (IH) and quasi- 
degenerate (QD) neutrino mass spectra one has 

HI — > mi <C m 2 , m 3 , Amg = AmJi ; 
IH — > mi <C m-2 , 777,3 ) Am?) = Am 2 2 ; 
QD — > mi ~ 7772 — 7773 , Anig = Am|i , (I ) 

where Am?, and Am 2 tm = Am^ ^> Am?, are the solar 
and atmospheric neutrino mass squared differences, re- 
spectively. At the la level, the allowed ranges for Am| 
and Am 2 tm are H3, lir| 

Am| = (6.5 - 8.5) x 1CT 5 eV 2 , 

Am 2 tm = (2.6±0.4)x I(T 3 eV 2 , (2) 

with the best-fit values 

Am% = 7.13 x ICT 5 eV 2 , 

Am 2 tm = 2.6x I(T 3 eV 2 . (3) 

For three light Majorana neutrinos, the leptonic mixing 
matrix U can be written as 

/ u el u e2 u e3 \ 

U=\ LV U^ 2 E/> =C/ <5 diag(e- M ,e^,l), (4) 

V u Tl u T2 u t3 J 

where a and (3 are Majorana-type phases and Us can be 
parametrized in the form 

/ C3C2 S3C2 s 2 e~ lS \ 

Us = -S3C1 - c 3 S!S 2 e lS C3C1 - s 3 s 1 s 2 e 1 ' s sic 2 

\ S3S1 - C3CiS2e i<5 -C3S1 - S3Cis 2 e 115 cic 2 / 

(5) 

with Si = sin^i, c; = cos^ (7 = I, 2, 3) and 6 the Dirac- 
type CP-violating phase. For Dirac neutrinos U reduces 
to Us, due to the absence of Majorana phases. Depending 
on the type of neutrino mass spectrum, the solar, atmo- 
spheric and CHOOZ [43 mixing angles (#©, 6> atm and 9, 



3 



respectively) can be extracted from U in the following 
way: for the HI and QD neutrino mass spectra, 



tan 6*0 = 
sin 6 



\U e ; 



= \U e3 \ = sin 6*2 



tan #3 , tan a 



\U, 



u. 



t:>,\ 



and for an IH spectrum, 

1^3 1 



tan0 m = 



tan( 



\U, 



ft i 



Prl\ 



sm f 



= tan 9\ , 
(6) 

Wel\. (7) 



In this case the expression of 0©, a t m and in terms 
of 6*i_3 m the parametrization of Eq. (jSJ) is not simple. 
From the global analyses performed in Refs. |4jj and |4t| . 
the solar and atmospheric mixing angles are constrained 
to lay in the la intervals 



tan 2 = (0.33 - 0.47) 
with the best-fit values 

tan 2 00 = 0.39 , 



sin 2 20 atm = 1.0018;$ . ( 8 ) 



sin 2 20 atm = 1.00. 



(9) 



electroweak scale, being A a 3 x 3 symmetric matrix of 
dimensionless couplings of order unity and A the scale at 
which this interaction is generated. These vertices can be 
depicted by the Feynman diagrams in Fig. ^ where the 
flavor dependence of each vertex is explicitly shown. The 
arrows in the scalar lines indicate the flow of the positive 
charge for the T — 1/2 component of the doublet. 







For sin0 we quote the result from combined analysis of 
the solar neutrino, CHOOZ and KamLAND data per- 
formed in Ref. 481, 



sin0 < 0.18, 



(10) 



with a 95% confidence level (CL). 

The upcoming long-baseline neutrino oscillation exper- 
iments will face the challenge of detecting CP- violating 
effects induced by the Dirac phase S The difference 
of the CP-conjugated neutrino oscillation probabilities 
P(v e — ► v^) — P(& e — * Vfx) is proportional to the quantity 



— - sin 20i sin 202 sin 203 sin 6 . 



(11) 



Present estimates indicate that for |J7"| > 10 -2 it will be 
possible to observe CP violation effects in these experi- 
ments. 



III. PARAMETRIZING THE NEW PHYSICS 
CONTRIBUTIONS TO LEPTON MASSES 



FIG. 1: Feynman diagrams corresponding to the vertices in 
Eq. Ill 2t . In each vertex the flavor dependence is explicitly 
shown. 

There are different mechanisms that may lead to the 
neutrino mass operator given in Eq. I|12|l . In particular, 
it may be originated by a heavy scalar triplet A, which 
is conveniently written in matrix form as 



A 



A+ V2A+ 4 
V2A° -A+ 



(13) 



This scalar triplet couples to the lepton and Higgs dou- 
blets as 



C A =Y i fl i &U 



^Pj + g^tfAiazf + H.c. , (14) 



where Y A is a matrix of Yukawa couplings and has 
dimension of mass. The exchange of the heavy triplet 
results in an effective neutrino mass operator, given by 



A 



= Y s 



A.9A0 



(15) 



Before electroweak symmetry breaking (EWSB), the 
terms of the Lagrangian that originate the charged lepton 
and light Majorana neutrino masses can be written as 



ia 2 I 



H.c. 



C = -Y£ l u e Rj + ^ {£ Li io 2 (j>*) 

(12) 

where l Li = (y Li e Li ) T , <j) = (0+ 4>°) T is the SM Higgs 
doublet and Y e is the usual 3x3 matrix of the charged 
lepton Yukawa couplings. The second term is a neu- 
trino mass operator |50j generated by physics above the 



with 771a the mass of the triplet. 

Another possibility is the exchange of heavy right- 
handed neutrinos vr (the seesaw mechanism). In this 
case, the relevant terms are 



r — 



1 ij *-Lz 



ia 2 <j>* VRj - ^M Rij v c Ri v Rj + H.c. , (16) 



with Y v is the Dirac neutrino Yukawa coupling matrix 
and M R the heavy right-handed neutrino mass matrix. 
The integration of the heavy neutrino fields generates an 
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effective mass operator for the light neutrinos of the form 



(17) 



After EWSB, the terms in Eq. itl"2l yield the mass 
terms for the charged leptons and left-handed neutrinos. 
Using matrix notation in flavor space, the mass terms 
can be written as 



C, 



-e L M e e R - -v L M L v c L , 



(18) 



where M e = vY e and M L = 2A v 2 /A, with v = 174 GeV 
the vacuum expectation value (VEV) of the SM Higgs 
boson. 

Our parametrization of the new physics effects is ob- 
tained after general considerations on weak basis inde- 
pendence. For convenience, let us define 

Ml=sA = ^M l = ^, (19) 

with e a dimensionless parameter to be specified later. 
Under the change of weak basis 



h = V L £' L , e R = V£e' R 



we have 



Y e -> Vl Y e V% , M L -> V[M L Vl 



tji> 



(20) 



(21) 



Let us assume that some perturbations arising from ra- 
diative corrections are added to the "tree-level" matrices 
Y e and Ml, 



Y' 



Y e + SY e 



Mr 



M L + SM L 



(22) 



The matrices SY e , 5 Ml are functions of Y e , Ml and 
other SM and new physics parameters. Under the change 
of basis defined in Eqs. I|2(JI) . the perturbations must 
transform as 



5Y e -> VI SY e V£ , 8M L -> VI SM L U L * 



(23) 



since the physical observables must be independent of the 
choice of weak basis. 2 These transformation laws imply 
that the perturbations have the form [2^ 



5Y e = X e Y e + ( e Y e Y e1 <Y e +r) e M L M L Y e + 
5M L = XlMl + ClMlMWl 



Y e Y e Hl L + M L Y e *Y eT ) +■■■ , (24) 



2 These transformation properties for 8Y e and 5Ml do not assume 
that the Lagrangian is invariant under the change of basis in 
Eqs. alone. Within the SM, the La grangian is invariant 

under these transformations, but this does not happen in some 
of its extensions, for instance in the MSSM. Besides, at very high 
energies, some flavor symmetry might single out a special weak 
basis. Below that scale, and in particular at low energies, this 
symmetry is broken. 



where the Aj, Q and Tji coefficients (i — e,L) are func- 
tions of Y e , Ml and of the coupling constants, which are 
invariant under the transformations of Eqs. (121 H and in 
general complex. The higher order terms in this expan- 
sion are expected to be smaller. The effect of the A, terms 
in Eqs. (12411 is to rescale the masses by common factors 
(1 + A e ) for charged leptons and (1 + Al) for neutrinos, 
without affecting neither the mass hierarchy nor the mix- 
ing. The Q terms also rescale the masses, but with a dif- 
ferent factor for each lepton: m ej — > m ej (1 + Ce rn2 ej / v ' 2 ) 
for charged leptons and rrij — > rrij (1 + (l m 2 /Af 2 ) for 
neutrinos. Hence, the Q terms modify the mass hier- 
archies. The rji terms are the lowest-order ones which 
modify the leptonic mixing. 

For some SM extensions, in principle there may exist a 
matrix X (not necessarily square) of couplings between 
the leptons and other particles, transforming as Y e or 
Ml either on the left or on the right side. For instance, 
if X transforms under the change of basis in Eqs. (|20J) as 



X 



(25) 



this matrix would contribute to 5Y e and 5Ml with terms 
XX^Y e and (XX* M L + M L X*X T ), respectively. Such 
possibility will not be considered in the following, and in 
this respect our analysis is not be the most general one. 3 
We thus assume Y e and Ml as the only sources of flavor 
violation in the lepton sector. This case can naturally 
arise if some symmetry relates the couplings in Eq. i|12|) 
with the ones between the leptons and the new particles. 

It is worthwhile showing some examples of Feynman 
diagrams which contribute to the different terms in 
Eqs. 124JI within the SM, including also the effective neu- 
trino mass operator in Eq. (|12f> as part of the SM vertices. 
The Xi terms result from diagrams with minimal flavor 
structure, as for example diagrams (a) and (b) in Fig. 
with the exchange of a B boson with flavor-universal cou- 
plings. The remaining terms in Eqs. (|24|l require the 
exchange of one or more <fi doublets. In particular, the 
Q e and t]l terms arise from diagrams like (c) and (d) 
in Fig. respectively. At the one loop level the terms 
with rj e and Cl are absent, and to generate them it is 
necessary to consider two-loop corrections, for instance 
diagrams (e) and (f), respectively. 

In new physics scenarios there are additional interac- 
tions that may or may not be suppressed by a large scale 
A'. These interactions mediate Feynman diagrams giving 
further corrections to the Il<\>&r and iL^L ( t ) * ( t'* vertices. 
Several examples of new physics contributions to these 
operators can be found in Ref. We remark that, if 

the particle content and parameters of the new physics 



3 Due to the ignorance regarding the structure of this matrix X, 
the discussion of the effects of the terms XX ( Y e , {XX^M L + 
M L X*X T ) is not feasible. 
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<l> 




(e) (f) 

FIG. 2: Examples of SM diagrams giving corrections to the 
Il<\>&r and £l^l4 > *4 > * vertices. In each vertex the flavor de- 
pendence is explicitly shown. 



model are specified, the corrections to the charged lep- 
ton and neutrino mass matrices can be completely de- 
termined. However, in the absence of any experimental 
indication favoring any of the theories beyond the SM, 
the Aj, Ci) Vi coefficients cannot be predicted. Then, it is 
sensible to perform a statistical analysis in order to de- 
termine, under certain assumptions, how much the tree- 
level predictions can change due to radiative corrections 
from new physics, parametrized according to Eq. (|24|l . 
leaving A^, r\i as unknown parameters. The size of 
these coefficients is expected to be similar, because all 
the terms in Eqs. H24J1 can be generated at the one loop 
level (although the terms with r/ e and C,l appear at next- 
to-leading order within the SM, they may arise at leading 
order in other models, e.g. with a scalar triplet). On the 
other hand, the higher order terms omitted in Eqs. Q24JI. 
involving products of 5 matrices or more, are expected 
to be suppressed by a factor ~ 10 with respect to the 
leading ones. 

One crucial issue for our analysis is the value of the 
parameter e in Eq. I|19|l . which accounts for the normal- 
ization of Ml- The size of this parameter reflects the 



suppression of the new interactions, and determines the 
relative importance of Ml with respect to Y e in the ex- 
pressions used for the perturbations. We consider two 
limiting scenarios: 

(i) In the first scenario we take e = 1, in which case 
Ml = A with matrix elements of order unity. This 
corresponds to a situation where there are new interac- 
tions which are not suppressed by a large scale A'. In 
this scenario, the terms C, L M L M\M L and r) e M L M L Y e 
in Eqs. I|24[) are not negligible, and have an important 
influence on the neutrino mass hierarchy and mixing, re- 
spectively. 

(ii) In the second scenario we assume e <C 1, so that 
these two terms (which have two or more powers of Ml) 
can be omitted in Eqs. J23J), resulting in 

6Y e \ S2 = A e Y e + C e Y e Y e ^Y e + ■ ■ • , 
SM L \s2 = Ai M L + vl (y c Y^ M L + M L Y e *Y eT ^j 

+ ■••. (26) 

This scenario corresponds to a new physics model in 
which the new interactions are suppressed by a large scale 
A' ~ A, for instance in models based in the seesaw mech- 
anism. In the limit e -C 1, the normalization of Ml is 
irrelevant. 

From Eqs. I|26|) we see that in scenario 2 the expres- 
sions for the perturbations are formally similar to the 
one loop RG equations within the SM. Therefore, in this 
scenario the results obtained within our framework are 
expected to be similar to the results of RG evolution, 
bearing in mind that in the case of the RG equations the 
coefficients Ai, Ce and r\L are fixed, while in our case they 
are unknown in principle. 

If the neutrinos are not Majorana but Dirac parti- 
cles, the Yukawa term of the Lagrangian originating their 
masses reads 

C D = -Y% hi ^C720* vrj + H.c. (27) 

For this term and the charged lepton one, the change of 
basis analogous to Eqs. (|20l) reads 

e L = V L l' L , e R = V£e R , vr = V^v' r . (28) 

Under this transformation, the Yukawa matrices trans- 
form as 

Y e -> Vl Y e Y v -> Vl Y v V£ . (29) 

This allows to obtain the expressions for the perturba- 
tions for Dirac neutrinos, 

SY e = \ e Y e + Q e Y e Y e ^Y e + r, e Y v Y v ^Y e + ■ ■ ■ , 
SY V = \ y Y v + Q u y v Y v * s Y u + r} v Y e Y e * s Y lJ 

+ ■■■. (30) 

In the following, we will generally refer to the case where 
neutrinos are Majorana particles, and quote the results 
for Dirac neutrinos when relevant. 
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IV. NUMERICAL RESULTS 

Using the parametrization of new physics contributions 
given in Eqs. (|24[l . Q26p. (|30|) . we study the changes in 
the pattern of neutrino masses and mixings due to these 
corrections. For this purpose, we take the unknown coef- 
ficients Xi, Ch r\i as random complex parameters, gener- 
ated with a Gaussian distribution centered at zero and, 
for simplicity, we assume that the standard deviations 
coincide: 



(|Ai 



(10 



K . 



(31) 



Contrarily to what could be expected, this is not a serious 
bias in the analysis, because the moduli of the random 
parameters are not constrained to be all equal, and only 
the standard deviations of the distributions are assumed 
to be the same. The phases of A^, Q and r]i are generated 
uniformly between and In. In our study the following 
procedure is applied: we fix a value of k and generate 
a large set of matrices using Eqs. <j2"4T) . it2"6*|) or i|3U|) . as 
appropriate, with random coefficients A^, £i: Vi- These 
matrices are diagonalized in order to obtain the masses 
and the neutrino mixing matrix. We then select some 
observable, and examine its distribution over the set of 
matrices. For each value of /t, the la limits on this ob- 
servable are defined as the boundaries of the 68.3% confi- 
dence level central interval, evaluated from the sample of 
random matrices. These la limits reflect the "average" 
behavior of the observable under consideration, when ar- 
bitrary perturbations are added to the original matrices. 
It must be emphasized that the maximum and minimum 
values can be very different from the average values, and 
some situations are found where in average the observable 
does not change appreciably under perturbations, but for 
fine-tuned values of the random parameters it does. 

In the numerical analysis we take initial "tree-level" 
matrices Y e and Ml (Y v for Dirac neutrinos) reproduc- 
ing the current experimental data summarized in Section 
ITT1 The charged lepton masses are taken at the scale Mz 
|52|. We assume sin6> = 0.15 and fix the Dirac and Ma- 
jorana phases to be S — 7r/2, a — 7r/3 and (3 = 7r/5. 
We analyze separately the two possibilities of a HI or a 
QD spectrum. For the case of an inverted hierarchy the 
results turn out to be very similar to those found for a 
normal hierarchy, and we do not present them. We sum- 
marize our input values in Table [H In scenario 1, we take 
M l with the mass of the heaviest neutrino normalized to 
unity. The normalization of Ml is irrelevant in scenario 
2, as shown in the previous section. 



A. Stability of mass and mixing parameters 

Let us discuss how the parameters r = Am^/Am^, 
tan 0q, tan# atm and sin# change when corrections are 
added to the mass matrices. Our aim is to investigate 
whether these quantities are stable or not, and to what 



Parameter 


Value 


m e 


0.487 MeV 


TUu 


0.103 GeV 


m T 


1.747 GeV 


riii 


J 10 -5 eV (HI) 


1 1 eV (QD) 


ArriQ 


7 13 x 10~ 5 eV 


Am 2 , 


2 6x1 -3 eV 


tan 00 


0.62 


T, r\ ,ff r7 n +- m 
uuiil ** ax in 


1 


sin 6 


0.15 


8 


tt/2 


a 


7r/3 


P 


tt/5 



TABLE I: Input parameters used for the unperturbed mass 
matrices. 



extent they get modified by the perturbations. We first 
present the results for scenario 1, and later discuss the 
differences with scenario 2 and the results for Dirac neu- 
trinos. 

Regarding the ratio of mass squared differences r, we 
observe in Fig.[3]that the corrections to the matrices have 
a large impact on this quantity, both in the cases of a HI 
or QD spectrum. This plot (and the remaining ones in 
this section) must be interpreted with caution: it does 
not provide any limit on the size of the corrections, on 
the basis of the experimental measurement of r, because 
the initial "tree-level" value we use for r needs not be 
equal to the observed value. Instead, the meaning of the 
plot is that r is not stable under perturbations, and from 
an initial value chosen to be r = 0.027 one can obtain 
values between 0.021 and 0.034, for a HI spectrum and 
K = 0.2. 




FIG. 3: Effect of the perturbations on the ratio r in scenario 
1. For illustration, the present la limits are displayed on the 
vertical axis. 

The effect of the perturbations on tan 9q is quite differ- 
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ent for a HI or QD spectrum. In the former case, tan 9q 
is very stable even for relatively large perturbations, as 
can be noticed in Fig. On the contrary, for quasi- 
degenerate neutrinos, the value of tan# Q can change sig- 
nificantly with new physics corrections. This fact sug- 
gests that, if neutrinos are quasi-degenerate, the under- 
lying tree- level pattern of lepton mass matrices could cor- 
respond to bi-maximal mixing, being the observed value 
tan 00 ~ 0.6 the result of radiative corrections. We will 
analyze in detail this possibility at the end of this section. 
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FIG. 4: Effect of the perturbations on tan 9q in scenario 1. 
For illustration, the present la limits are displayed on the 
vertical axis. 

The behavior of tan a t ro is the opposite to the one ob- 
served for tan 0q, as it can be perceived from Fig. [3] for a 
HI spectrum this parameter is modified by perturbations 
on the mass matrices, while for a QD spectrum it is fairly 
stable. This shows that, for the case of quasi-degenerate 
neutrinos, the experimental observation of tan0 atm ~ 1 
must correspond to tan a tm — 1 in the mass matrices, 
because this prediction is not altered by the corrections. 
On the other hand, for a HI spectrum the observation of 
tan f? atm ~ 1 could be either a coincidence, or result from 
a specific symmetry leading naturally to this value and 
making higher order corrections very small. 

The analysis of sin 9 (which equals | U e 3 1 for a normal 
hierarchy) shows that it does not change under perturba- 
tions for a QD spectrum, but it is considerably modified 
when the neutrino masses are hierarchical (see Fig. IBfr. 
However, one striking feature of our analysis is that the 
ratio 



sin 9 
tan 6> at , 



(32) 



remains practically constant even when large perturba- 
tions are added to Y e and Ml- for k = 1 sin (9 and 
tanf? atm change by more than ±40%, while their ratio 
changes less than 1%, as can be noticed in Fig. 

This feature may have important consequences for 
model building. When sin0 is experimentally measured, 
the ratio ui will be an excellent tool to investigate the 




FIG. 5: Effect of the perturbations on tan# a tm in scenario 
1. For illustration, the present la limits are displayed on the 
vertical axis. 
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FIG. 6: Effect of the perturbations on sin# in scenario 1. 



0.17 r 



0.16 - 



3 0.15 - 



0.14 



13 1 1 1 L 

0.2 0.' 



Hierarchical 
Quasi-degenerate 



0.6 0.8 



FIG. 7: Effect of the perturbations on the ratio u 
sin#/tan# a t m in scenario 1. 



structure of the lepton mass matrices, because it is very 
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insensitive to radiative corrections (as long as these cor- 
rections do not have any additional source of flavor vio- 
lation, which is the framework we consider). This ratio 
will then allow to test experimentally, with a "clean" ob- 
servable, the textures for lepton mass matrices implied 
by flavor symmetries proposed in the literature. 

The remaining parameters to be investigated are the 
CP violating phases, for which the results depend on 
the initial values used (see Table QJ. For a HI spectrum, 
the Dirac phase S remains virtually constant at its initial 
value, while the Majorana phases a and j3 vary over a 
wide range, 0.6 < a < 2, 0.4 < (3 < 2.2 for n = 1. 
For a QD spectrum, the Majorana phases remain within 
±10% of their initial value for k = 1, while the Dirac 
phase varies in the interval 1 < 5 < 1.8. 

In scenario 2, the behavior is very different for a HI 
spectrum. In this case, we find that r, tan0 Q , tan# atm , 
sm8 and the three CP violating phases remain constant 
when the perturbations in Eqs. (j26|) are added to the 
matrices. In the SM these quantities exhibit a similar 
behavior under RG evolution |53[ . as its expressions are 
formally identical to ours. On the other hand, in the case 
of a QD spectrum, the differences between scenarios 1 
and 2 are not significant, and the discussion above applies 
also to scenario 2. This contrast can be understood in 
view of the analysis of the dependence on the neutrino 
masses presented in Section HV Bl below. 

For Dirac neutrinos, the results are found to be rather 
similar to the ones obtained in scenario 2. For a HI spec- 
trum the parameters under consideration do not change 
when the perturbations in Eqs. (|30|l are included. For 
a QD spectrum, however, they are modified, following 
practically the same behavior that can be observed in 
the plots for scenario 1 shown in this section. 

Finally, we note that the results on r, tan(9 Q and 
tan# a tm are almost independent on the value of sin# 
used. Other quantities obviously depend on the particu- 
lar value of sin#, as for example J , which is proportional 
to sin 9. The choice of the CP- violating phases is only rel- 
evant for the "tree- level" values of quantities that depend 
on them, like J and m ee . For a QD spectrum, the devi- 
ations of r, tan 00, tan# atm , sin# and u> are very similar, 
and for a HI spectrum the influence of phases on these 
quantities is completely negligible. We have also checked 
that our results do not change when the next terms in the 
expansion of Eqs. I|24|) (with products of 5 matrices) are 
included, even in the unrealistic limit where these terms 
have similar coefficients. We have found that the quanti- 
ties that are stable remain stable, and the quantities that 
change under perturbations exhibit an analogous behav- 
ior with the inclusion of these terms. 



B. Dependence on the neutrino masses 

We have verified that the influence of the perturba- 
tions on some parameters depends strongly on the type 
of neutrino spectrum, namely the deviations of tan# Q 



are negligible for a HI spectrum while they are large if 
the neutrinos are quasi-degenerate. It is then convenient 
to analyze the dependence of the deviations on the mass 
of the lightest neutrino, which we take between 10 -5 eV 
and 2.2 eV (the direct bound from the Mainz and Troitsk 
experiments). In Figs.lSHTTlwe plot the effect of the cor- 
rections on r, tan 00, tan# atm and sin#, respectively, for 
k = 0.2. 

0-08 1 




m, (eV) 



FIG. 8: Effect of the perturbations on r for k = 0.2, as a 
function of the mass of the lightest neutrino, in scenario 1. 




m, (eV) 



FIG. 9: Effect of the perturbations on tan 8q for k — 0.2, as 
a function of the mass of the lightest neutrino, in scenario 1. 

It is apparent that these parameters are remarkably 
more stable in the region around mi = 0.3 eV than 
for the rest of values of mi. This can be understood 
as follows: the only two terms that influence the mix- 
ing are Tj e M L M f L Y e and rj L {Y e Y e ^M L + M L Y e *Y eT ^j . 

Of these, the former is relevant only for a HI spectrum, 
whereas for a QD spectrum it does not have any influ- 
ence. On the contrary, the latter term is important for 
a QD spectrum but its impact is negligible if the neu- 
trino masses are hierarchical. Hence, the deviations in 
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FIG. 10: Effect of the perturbations on tan# atln for k — 0.2, 
as a function of the mass of the lightest neutrino, in scenario 
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FIG. 12: Effect of the perturbations on uo for k = 0.2, as a 
function of the mass of the lightest neutrino, in scenario 1. 
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FIG. 11: Effect of the perturbations on sin# for k = 0.2, as a 
function of the mass of the lightest neutrino, in scenario 1. 



tan^atm and sin (9 appearing in the left part of Figs. I1UI 
and 1111 are due to the r\ e term, and the deviations in R, 
tan 0q, tan0 a t m and sin# that can be seen in the right 
side of Figs. ISHTTl are a consequence of the t\l term. The 
deviations in r in the left hand side of Fig. [SJ is an effect 
of the Ci term, which does not contribute to the mixing. 

The region m ~ 0.3 eV is of special interest, since these 
neutrino masses will be probed in forthcoming experi- 
ments like KATRIN, which is planned to start in 2007. 
If the mass of the lightest neutrino happens to be in this 
range, this will mean that the corrections to the tree-level 
mass matrices will have a much smaller impact on the hi- 
erarchy of mass squared differences and the mixing. The 
same is also true for the Dirac and Majorana phases. For 
completeness, in Fig. Elwe show the la limits on u> for 
the same range of mi . We observe that this ratio remains 
virtually constant in the whole interval. 

For scenario 2 and for Dirac neutrinos, the dependence 



on the lightest neutrino mass is much simpler. For a HI 
spectrum, all the quantities studied are stable, and in 
these cases the effects of perturbations are as in Figs. 00- 
^]but without the deviations present in the left part of 
some of these plots. 



C. Stability of flavor relations 

We are interested in finding flavor relations which are 
stable under perturbations of the mass matrices. By "sta- 
bility" we mean that, if these relations hold for the tree- 
level matrices, they still hold to a good approximation 
when perturbations are added. With some exceptions, 
most flavor relations found in the literature correspond 
to Majorana neutrinos with a HI spectrum. In scenario 
2, the parameters r, tan0 Q , tan# atm and sin 5, as well 
as the CP violating phases, remain constant for a HI 
spectrum. Therefore, any flavor relation among these 
parameters is stable. In scenario 1, most of the flavor re- 
lations studied are affected by the random perturbations. 
For illustration, we show the effect of perturbations on 
some relations. For texture Al in Ref. |54|. there are two 
predictions, 



1 



sin 

m 



TOatm 



= -tan0 atm sin20 \/r, 



sin B a Vr, 



(33) 



where m ee is the effective mass for n eutrinol ess double 
beta decay processes, and m a t m = \J Am atm . In order 
to test the stability of these relations under the radiative 
corrections considered here, we define the ratios 



Ri 
Ru 



1 tan 9 a 



■ sin 29, 



sin 



sm( 
2 a 



\m ee \/m a .tn 



(34) 
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which equal unity for the tree-level matrices. These ratios 
are plotted in Figs. ^| and ^] respectively, where we 
have used sin# = 0.074, S = a — 0, [3 = 0.85, so that 
the initial matrices fulfill these relations. The rest of 
parameters is taken from Tabled 




FIG. 13: Effect of the perturbations on the ratio Ri, in sce- 
nario 1. 




FIG. 14: Effect of the perturbations on the ratio Ru, in sce- 
nario 1. 

From Fig. ^| we see that the first relation is modified 
when corrections are added to the mass matrices. At 
any rate, the deviations on this relation are much smaller 
than the changes in r, tan^atm and sin# (see Figs. |31 
EJ). We also notice from Fig. ^] that the accuracy of the 
second relation is hardly affected by perturbations on the 
mass matrices. This feature makes its experimental test 
cleaner, and less dependent on unknown corrections to 
the tree-level textures. For texture Bl of Ref. [54| (see 
also Ref. £3) we have 



1 



sin l 



= — tan 6> atm tan 20© \J r cos 20(. 



(35) 



which differs from the first of Eqs. (|33[1 by factors de- 
pending on 9q . Since 9q is stable for a HI spectrum, the 



effect of corrections on this relation is similar, and the 
plot obtained is identical to Fig. ^| Another interesting 
relation is 56] 



sin ( 



(36) 



The left-hand side of this equation varies with the per- 
turbations, but the right-hand side does not. We define 



_ 1 frn^ 
mu - -r—z\ — 
sin u y m M 



(37) 



and set sm9 = 0.0688 (with the rest of parameters as in 
Table[IJ in order to test the stability of this relation. The 
result can be seen in Fig. 1151 




FIG. 15: Effect of the perturbations on the ratio Rui, in 
scenario 1. 

The conclusion one may draw from the study of these 
examples is the following: if the new physics interactions 
are suppressed (this situation corresponds to scenario 2), 
the flavor relations are stable and, if they hold at tree- 
level, they also hold when the corrections in Eqs. (|26|l are 
included. On the other hand, if the new interactions are 
not suppressed and the corrections have the full form 
of Eqs. (|24|l . with Ml of order unity (this possibility 
corresponds to scenario 1) some of these flavor relations 
are modified by the perturbations. 



D. Special limits 

We have previously remarked that, for quasi- 
degenerate neutrinos, the perturbations in the mass ma- 
trices modify tan#© but do not affect tan# atm signifi- 
cantly. Then, one interesting question naturally arises: 
Is it possible to have a QD spectrum with bi-maximal 
mixing at the tree level, so that the smaller value of 
tan 0q ~ 0.6 is due to effects of new physics? To test 
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this hypothesis, we set tan6* = tan# atm = 1 in our ma- 
trices, with the rest of parameters as in TableQJ and ana- 
lyze how tan 9q and tan # a t m change when perturbations 
are added. For scenario 1, the results are displayed in 
Fig. For scenario 2, the results are shown in Fig. IT7I 




FIG. 16: Effect of the perturbations on tan#0 and tan# atm , 
for quasi-degenerate neutrinos and initial bi-maximal mixing, 
in scenario 1. 




FIG. 17: Effect of the perturbations on tan#© and tan# a tm, 
for quasi-degenerate neutrinos and initial bi-maximal mixing, 
in scenario 2. 

From these figures we conclude that in both scenarios 
it is possible that, from an initial bi-maximal pattern, 
large corrections to the mass matrices modify signifi- 
cantly tan 0q, bringing it to its experimental value, while 
keeping tan# atm ~ 1. This result is only slightly depen- 
dent on the value of sin# used, and is valid for mi > 0.6 
eV. For Dirac neutrinos the same effect is found, and the 
results are very similar to the ones in scenario 2. 

Other interesting situation corresponds to sin 8 = at 
the tree level. In this case, for a QD spectrum a nonzero 
sin 6? can be generated by the perturbations. If the neu- 
trino masses are hierarchical, the value of s'm9 induced 



by the corrections is negligible. The results for scenarios 
1 and 2 are displayed in Figs.ll8land ll9l respectively. For 
Dirac neutrinos, the value of sin (9 generated is one order 
of magnitude smaller. 
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FIG. 18: Effect of the perturbations on sin#, for an initial 
vanishing value, in scenario 1. 
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FIG. 19: Effect of the perturbations on sin#, for an initial 
vanishing value, in scenario 2. 

Finally, we consider the situation when 5 = at the 
tree level but the Majorana phases are not zero. In this 
limit, the CP violating parameter J in Eq. Ijllfl van- 
ishes. For a QD spectrum, the corrections to the mass 
matrices induce a Dirac phase (provided at least one of 
the Majorana phases is nonzero) large enough to yield 
J r~j 10~ 2 , which may be observable by future long base- 
line neutrino oscillation experiments [4!j ■ This is shown 
in Figs. E0| and El and holds for mi > 0.7 eV. Setting 
one of the initial Majorana phases to zero does not elim- 
inate this effect, and for a — ir/3, (3 — the values of J 
obtained are up to 0.02, even larger than the ones shown 
in Figs. 1201 and 1211 For a HI spectrum, the phase S gen- 
erated is negligible. In the case of Dirac neutrinos, only 
the Dirac phase is physically meaningful, and the phases 
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of the random parameters , Q and rji are not enough 
to produce a relevant value of J . 
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FIG. 20: Effect of the perturbations on J , for an initial van- 
ishing Dirac phase, in scenario 1. 
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FIG. 21: Effect of the perturbations on J , for an initial van- 
ishing Dirac phase, in scenario 2. 



V. OUTLOOK 

In this paper we have studied the possible effect of un- 
known corrections from new physics on the neutrino mass 
hierarchy, mixing and CP violation at low energies. We 
have focused on the case when neutrinos are Majorana 
particles, but have discussed the results for Dirac neutri- 
nos as well. We have proposed a general parametrization 
of the corrections to the tree-level mass matrices, based 
on weak basis invariance. Using this parametrization, we 
have examined the consequences of adding random per- 
turbations to the mass matrices, as a means to explore 
the possible effects that radiative corrections might yield. 

We have analyzed the stability against corrections of 
the ratio of mass squared differences and the mixing 



angles, for a hierarchical or quasi-degenerate neutrino 
spectrum. We have found that these quantities are 
generally modified by the perturbations, but the ratio 
lu = sin 9/ sin # atm is remarkably stable even under large 
corrections. This desirable property makes this quantity 
specially suited for the experimental test of specific tex- 
tures of neutrino mass matrices, since it is hardly mod- 
ified by unknown corrections to the tree-level matrices. 
We have also examined the stability of some flavor rela- 
tions predicted by models in the literature. 

The dependence of the deviations on the neutrino spec- 
trum has also been investigated. We have found that the 
region of neutrino masses m ~ 0.3 eV is specially stable. 
For neutrino masses around this value, the possible devi- 
ations in r, tan 0q, tan# atm , sin^ and the CP violating 
phases are rather small. This mass region is of special 
interest, since it will be tested in upcoming experiments. 

We have addressed the question whether the tree-level 
mass matrices could correspond to bi-maximal mixing, 
being the observed value tan# Q ~ 0.6 the result of ra- 
diative corrections. We have demonstrated that, from 
an initial bi-maximal pattern, in the case of a QD spec- 
trum the corrections to the mass matrices can bring down 
tan 8q to its experimental value while keeping tan 6 aiul 
very close to unity. This result holds both for Majorana 
and Dirac neutrinos, and for a lightest neutrino with a 
mass larger than ~ 0.6 eV. Therefore, the possibility of 
bi-maximal mixing cannot be excluded for a QD spec- 
trum. On the other hand, for a HI spectrum tan# a t m 
is modified but tan# not, then bi-maximal mixing is 
highly unlikely in this case. 

Other interesting limit examined is when sin 8 = at 
the tree level. In this case, large corrections to the mass 
matrices could yield sin9 ~ 10~ 2 if the neutrinos are 
quasi-degenerate. On the contrary, we have shown that 
for a HI spectrum or Dirac neutrinos the value of sin# 
generated by perturbations is negligible. 

Finally, we have investigated the situation when the 
Dirac phase in the mixing matrix vanishes at the tree 
level. In this case, the Majorana phases present can in- 
duce a non-vanishing Dirac phase in the mixing matrix 
by means of the perturbations. This Dirac phase is large 
enough to yield J > 10~ 2 , leading to observable CP 
violation effects in long baseline neutrino oscillation ex- 
periments. 
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